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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

LOW-AMPLITUDE DAMPING-IN-PITCH CHARACTERISTICS OF FOUR 
TAILLESS SWEPT WING-BODY COMBINATIONS AT MACH 
NUMBERS FROM O .85 TO 1.50 AS OBTAINED 
WITH ROCKET-POWERED MODELS 
By Charles T. D'Aiutolo 

SUMMARY 


A free-flight investigation was conducted to determine the damping- 
in-pitch characteristics of four rocket-powered wing-hody combinations 
between the Mach nmhers of O .85 and I. 30 . All models were tested at the 
same center-of -gravity location (17 percent mean aerodynamic chord) and 
had identical bodies thereby allowing a direct comparison of the aerody- 
namic characteristics of the various wings. The wings tested were an 
aspect-ratio-4, 45° swept composite-plan-form wing (derived from an 
aspect-ratio-6 wing by increasing the area behind the trailing edge of 
the Inboard 40 percent of the wing), an aspect-ratio-4, 45° swept highly 
tapered wing, an aspect-ratio-5, 52 . 5 ° swept highly tapered wing, and an 
aspect-ratio-5, 60° swept highly tapered wing. The composite wing had 
NACA 65 AOO 6 airfoil sections parallel to the free stream, whereas the 
swept highly tapered wings had 65A004 airfoil sections parallel to the 
free stream. 

The results indicated that, although all models were statically stable 
longitudinally throughout the Mach number range investigated, the dynamic 
stability of the models was low and differed appreciably for each model. 

The total damping factor and the slope of the lift curve were about the 
same order of magnitude resulting in values of the rotational damping-ln- 
pitch derivatives that were either negative (stable) or positive 
(unstable) . 

A comparison of the damping was made between a delta wing in combina- 
tion with a body and a swept wing in combination with a body and the 
results indicated that at M = I.I 6 the delta-wing configiiration had better 
damping characteristics . 
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INTRODUCTION 


Recently a series of systematic experimental investigations has been 
conducted on tailless aircraft configurations to verify the theoretical 
predictions of low dynamic longitudinal stability and regions of dynamic 
instability of these designs (see, for example, refs. 1 to 5). Inasmuch 
as most of this work was concerned with the delta-wing configurations and 
very little information is available concerning swept-wing configurations, 
this investigation was undertaken to determine experimentally the regions 
of dynamic instability at transonic speeds and low supersonic speeds of 
some tailless swept-wing airplane conf igurations . 

This report is a continuation of the investigation reported in ref- 
erence 5 contains the results from the flight tests, conducted by the 
Pilotless Aircraft Research Division, of four rocket -powered wing-body 
combinations. Two of the wings used in these configurations had 45^ of 
sweep and aspect ratios of 4, the only difference being that one was of 
composite plan form (resulting from a search for an improved design suit- 
able for high-speed flight, see ref. 6 ) and the other was highly tapered. 
The third and fourth wings were both of aspect ratio 3 and highly tapered, 
with one having 52.5^ of* sweep and the other 6 qo. The composite wing had 
NACA 65 AOO 6 airfoil sections parallel to the free stream, whereas the 
swept highly tapered wing had NACA 65A004 airfoil sections parallel to 
the free stream. The data are presented over a Mach number range of about 
0.85 to 1.50 corresponding to a Reynolds number range (based on respective 
mean aerodynamic chords) of about 7 to l4 X I 06 , respectively. 

The models were flown at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. The static and dynamic longitudinal 
stability characteristics of the models were determined by analyzing 
the oscillations produced in pitch by firing small pulse rocket motors 
which were mounted to provide thrust normal to the longitudinal axis of 
the models. The drag characteristics of the models and the stability 
characteristics were determined from the deceleration portion of the 
flights. 


SYMBOLS 

V velocity of flight, ft/sec 

V 

M Mach number, 

Speed of sound 

S total wing area, sq ft 
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c 

A 

A' 

I 

X 

P 

R 

q. 

b 

P 

k 


Cl 


g 


Cc 


g 

w 

Cjn, 


mean aerodynamic chord, ft 
aspect ratio 

cross-sectional area of any longitudinal station of complete 
model, sq, ft 

body length, ft 

distance along body measured from nose, ft 
air density, slugs/cu ft 

Reynolds number (based on mean aerodynamic chord of respective 
models ) 

free-stream dynamic pressure, Ib/sq ft 

total damping factor (logarithmic decrement of pitch oscillations) 
period of short-period oscillations, sec 

reduced-frequency parameter (based on respective mean aerodynamic 

chords of models ) , — 

2V 

lift coefficient, Cjj cos a - Cq sin a 

nonnal-force coefficient, — 

g q 

acceleration normal to reference axis as obtained from 
accelerometer 

chord-force coefficient, - ^ 

g q 


acceleration along reference axis as obtained from accelerometer 
positive for\ira,rd 


model weight, lb 


slope of lift curve per degree. 



Sa 


pitching-moment coefficient, P . 4 . t . q h lng moment 

qSc 
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Sc 

static -stability derivative per degree, 


Cmq - 


Sc, 


m 


Hi) 


, per radian 


SCjji 

Cm. = — ; Tf per radian 


2Y J 

C'Q drag coefficient, Cq cos a 4* Cj^ sin a 

a angle of attack, measured from flight path to fuselage reference 

line, deg 


& = 


1 da 


57.5 dt 


— , radians/sec 


angle of pitch, measured from horizontal to fuselage reference 
line, radians 


6 = — , radians/ sec 
dt 

t time, sec 

a) frequency of short-period oscillation, radians/sec 

[1 = relative -density factor 

pSc 

ly moment of inertia in pitch, slug-ft^ 

m mass, slugs 


MODELS AND APPARATUS 
Model Description 


The general arrangements of the models are shown in figure 1, and 
a photograph of one of the models is presented in figure 2. Geometric 
characteristics of the models are presented in table I. Model 1 had an 
aspect-ratio-4, 45° swept composite plan-form wing (derived from an aspect- 
ratio-6 wing by. increasing the area behind the trailing edge of the inboard 
40 percent of the wing) with a taper ratio which varied from about O .7 for 
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the outboard panel to about 0.5 for the Inboard panel and incorporating 
NACA 65 AOO 6 airfoil sections parallel to the free stream. Model 2 had 
an aspect-ratio-4, 45^ swept highly tapered wing, model 5 had an aspect- 
ratio-3, 52 . 5 ^ swept highly tapered wing, and model 4 had an aspect- 
ratio-5, 60 ^ swept highly tapered wing. The wings of models 2, 5 , and 4 
incorporated NACA 65A004 airfoil sections parallel to the free stream. 

Note that the wings of models 1 and 2 were identical in aspect ratio and 
sweep, but differed in plan-form shape and section thickness, whereas the 
wings of models 5 and 4 were identical except for wing sweep. 

Each model consisted of a basic fuselage to which was attached the 
wing to be tested. The fuselage was essentially a body of revolution 
with rather large wing-root fairings for structural purposes and consisted 
of an ogival nose section which contained the telemeter and a cylindrical 
body section which contained the wing mount, necessary fairings, the 
vertical tail, and the sustainer rocket motor. Construction of the fuse- 
lage was of aluminum alloy with magnesium skin. 

The wings of model 1 were constructed of wood with sheet-steel inlays 
(for structural purposes), whereas the wings of models 2, 5, and 4 were 
constructed of solid alumimmi. All wings were mounted on the fuselage 
(as shown in fig. 1) with the resultant center of gravity located at 
17 percent of the mean aerodynamic chord during the deceleration portions 
of the flights. 

Each model contained a Cordite sustainer rocket motor and was boosted 
by a lightweight ^--Inch HVAR rocket motor. All models with their 
boosters were launched from a rail-type launcher (as shown in fig. 5 ) at 
an angle of approximately 45^. 

The cross -sectional-area distribution of each of the models is 
presented in figure 4, for possible correlation of drag results. 


Ins t r ument a t i on 

Each model contained a standard 4 -channel NACA telemeter which trans- 
mitted continuous flight measurements of angle of attack (measured by a 
vane-type instrument located on a sting forward of the nose of the 
models), normal acceleration of the center of gravity, longitudinal accel- 
eration, and pitot stagnation pressure (measured by a tube located on a 
strut below the fuselage of the models). (See fig. 1.) 

The position of the models in space was determined from an SCE 584 
tracking radar and the velocity of the models was obtained by use of the 
Doppler velocimeter radar. Atmospheric data were obtained from a radio- 
sonde released just prior to each of the model flights. Fixed and 
tracking motion-picture cameras were used to observe the conditions of 
the models during the major portion of the flights. 
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TEST AND ANALYSIS 
Test 


The data for each model were obtained during the decelerating portion 
of the flights. The damping- in-pitch data were determined by disturoing 
each of the models in pitch by a series of four small rocket motors pro- 
viding thrust normal to the longitudinal axis of the model and located 
near the rear of the model. These rocket motors were timed to fire in 
sequence during the decelerating portion of the flights; however, the 
first small rocket motor to be fired in each of the four models fired 
prematiirely during the latter part of sustainer motor burning. Also, 
one of the remaining three small rocket motors in model L failed to fire, 
so that damping data from the flight test of model L were obtained from 
the firing of only two of these small rocket motors. Thus, in most 
instances, reliable damping data co\ild be only obtained from three of 
the disturbances. 

Time histories of angle of attack, lift coefficient, and Mach number 
covering the decelerating portions of the flight and the times at which 
the small pulse rockets were fired for each of the four models are shown 
in figure 5* Also shown are the envelopes drawn for each of the oscil- 
lations in pitch that were caused by the firing of the small rocket motors. 
The static and dynamic longitudinal stability derivatives were obtained 
from these oscillations. 

Data, obtained from spinsonde records, indicated that for each of 
the four models the rate of roll was approximately zero. 

Figures 6, 7> and 8 present, respectively, the variation of air 
density, velocity, and dynamic pressure with Mach mmiber for each of the 
tests. These quantities are presented so that a possible correlation of 
the data obtained from these tests with data obtained from other tests 
may be made. The range of Reynolds numbers of the tests is plotted against 
Mach number in figure 9 where the Reynolds number is based on the respec- 
tive mean aerodynamic chords of the models. 


Acciiracy 

The max i mum probable errors in the basic coefficients and angle of 
attack are as follows: 

M = 0.85 M = 1.25 


ACn 0.007 0.005 

ACc 0.007 0.005 

Aa, deg 0.20 0.20 
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These errors are largely due to errors in instrumentation which is a 
fimction of the full-scale range of the instruments and errors due to 
flight conditions. The Mach number as determined from Doppler radar 
velocity and radiosonde measurements is accurate to better than 2 percent. 


The analysis of the motions of the models presented in this report 
are based on the assumption of constant coefficients in the differential 
equations of motion. (See ref. 7 .) 

The slope of the lift curve CLq^ was found by plotting Cl against 
a as obtained from the flight records during an oscillation and graphi- 
cally measuring the slope. The Cm<x data that are presented in this 
paper were calculated by use of the expression 


which is based on a single-degree-of -freedom analysis instead of the more 
usual expression 


which is based on a two-degree-of -freedom analysis. For tailless con- 
figurations, the contribution of the total damping factor b to Chiq^ 
is negligible when compared to the frequency contribution, so that the 
single-degree-of -freedom expression given (eq. (l)) allows an accurate 
determination of Cm^. Values of the total damping factor b were 
determined by the method presented in the *’Analysis" section of refer- 
ence 5 ^ while values of Cmq + were determined by the method pre- 

sented in appendix A of reference 7 . 

The maximum value of the reduced-frequency parameter k = ^ (based 

on respective mean aerodynamic chords) that was determined from these 
tests was 0.025. Since this value is small, it is believed that the 
effects of the frequency of the oscillations of the models on the damping 
in pitch are not important in the determination of the damping-in-pitch 
derivatives and that the method of reducing the data used in this paper 
gives good results for these derivatives. IMpublished calculations, 
comparing the damping in pitch computed for terms to the order of k5 
with the damping in pitch computed for terms to the order of k, indicate 
that for k < 0.025 the difference is about 1 percent. These results 


Analysis 



(1) 



( 2 ) 
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indicated that terms of higher-order frequency are not required in the 
estimation of the damping-stability derivatives. 


RESULTS AND DISCUSSION 


The stability parameters of the models presented in this paper were 
determined from the coasting phase of the flights. Each of the models 
was tested with the center of gravity located at I7 percent of the mean 
aerodynamic chord behind the leading edge of the mean aerodynamic chord. 


Trim 

The trim characteristics of each model are shown in figure 5* 

Model 5 (a = 5, 52.5° sweep) experienced a small trim change in the tran- 
sonic speed range; however, it was the largest trim experienced by any of 
the models tested. At high subsonic speeds, model 5 flew at a trim eingle 
of attack of approximately 1.0° corresponding to a trim lift coefficient 
of approximately 0 .02 and experienced a nose-down trim change at transonic 
speeds. At supersonic speeds, model 3 flew at a trim angle of attack 
of approximately 0.5° corresponding to a zero-lift condition. 


Lift 

The variation of the slope of the lift curve for the four models is 
presented in figure 10. The test points represent data that were deter- 
mined from the oscillations produced by the firing of the small pulse 
rockets and, inasmuch as limited data were obtained, no attempt was made 
to fair curves through the test data. 

Models 1 and 2 (identical aspect ratio and sweep, but different in 
plan-form shape and section thickness) show about the same variation with 
Mach nimiber; althoiigh, at transonic speeds and supersonic speeds, CLa 
of model 2 is somewhat lower than the CLa °^ model 1. These differences 
are not due to wing flexibility since preflight static wing twist indi- 
cated that models 1 and 2 had about the same degree of flexibility at 
transonic and supersonic speeds. 

The effect of sweep on CLa ™s.y be seen by comparing model 5 
(52.5° sweep) with model 4 (60° sweep) since, except for wing sweepback, 
the models were the same. Also, both models had about the same degree 
of wing flexibility. The CLa of model 4 are somewhat lower than 

the CLa data of model 3 as would be expected since the wing sweep of 
model 4 was greater than the wing sweep of model 3- 
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Static Longitudinal Stability 

The static longitudinal stability for each of the four models is 
presented in figure 11 where it is seen that Cm^ increases with 
increasing Mach number through the transonic speed range, then decreases 
somewhat as the Mach number becomes supersonic. From a comparison of 
models 1 and 2 (same aspect ratio and sweep, but different in plan-form 
shape and section thickness), it is seen that at high subsonic and tran- 
sonic speeds model 1 has greater static longitudinal stability than 
model 2, while at supersonic speeds the reverse is true. A comparison of 
model 3 vith model k (same aspect ratio, but different sweepback) shows 
that model 5 bas slightly greater static longitudinal stability at high 
subsonic and transonic speeds than model 4, but at supersonic speeds the 
reverse is true. Note that the static stability of models 2 and 4 
increase to a higher Mach number before decreasing than models 1 and 5- 

Since all wings were mounted on identical bodies and were tested at 
the same center -of -gravity location with respect to the mean aerodynamic 
chord, a comparison of the aerodynamic-center location of the various 
wings can be made. These values of the aerodynamic -center location appear 
in the following table: 


Mach number 

Aerodynamic -center location in percent 
mean aerodynamic chord for - 

Model 1 

Model 2 

Model 5 

Model 4 

Subsonic .... 

0.361 

0.312 

0.341 

0.532 

Transonic . . . 

.564 

.424 

.325 


Supersonic . . . 

.415 

.464 

.492 

.492 


Note that generally for all models as the Mach number increases the 
aerodynamic -center location moves rearward; however, model 3 indicates 
a slight forward movement as the Mach number increases from high subsonic 
speeds to transonic speeds. It is noted that this forward movement of 
the aerodynamic -center location occurred while model 3 experienced a 
nose-down trim change (see fig. 5(c)). 


Dynamic Longitudinal Stability 

Total damping factor .- The total damping factor b was obtained 
from the time histories of the pitch oscillations (see fig. 5). This 
total damping factor includes the contribution of moment due to motion 
along a curved path at constant angle of attack 0^^, the moment due to 
plunging motion (vertical acceleration) and the translational 
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effect of CLq^. The variation of b with Mach number for each of the 
four models is presented in figure 12 . By use of the method presented 
in the "Analysis" section of reference 5 it was possible to determine 
numerous instantaneous values of b at high subsonic speeds for model 1 
and at low supersonic speeds for models 2, 3, and such that the solid 
curves presented in figure 12 represent values that were determined from 
experimental test data points only. The dashed curves represent fairings 
between other data points that were obtained for models 1, 2, and 3. 

Figure 12 shows that plan-form shape has a pronounced effect on the 
total damping factor for the low amplitudes reported in this paper. A 
comparison of models 1 and 2 (identical aspect ratio and sweep, but dif- 
ferent in plan-form shape and section thickness) shows that for the 
composite-plan-form configuration (model 1) the total damping factor b 
decreases abruptly then increases abruptly for small increments of Mach 
number at high subsonic speeds, while for model 2 this abrupt decrease 
followed by an abrupt increase in b over a small range of Mach nijmbers 
occurs at low supersonic speeds. These abrupt changes were determined 
from amplitudes that were less than 0.5° in angle of attack and therefore 
may not be significant. 

The total-damp Ing-f actor data for model 3 ( 52 - 5 ° sweep) and model 4 
(60° sweep) are also shown in figure 12. Inasmuch as limited data were 
obtained, it was not possible to determine the effect of sweep on the 
total damping factor. 

Rotational damping- in-pitch derivatives . - The rotational damping-in- 
pitch derivative for each of the foinr models was determined by the fol- 
lowing expression: 


where m is the mass of the model in slugs. It is seen from the fore- 
going equation ( 3 ) that the determination of the value of the rotational 
damping- in -pitch derivatives is dependent upon differences in the total 
damping factor b and the translational effect due to C]^. In the model 

tests reported in this paper (tailless configurations) the total damping 
factor and the translational effect of CLq^ are about the same order of 
magnitude resulting in small values of Cmq + Cm^ that are either nega- 
tive (stable) or positive (unstable). (See figure 13 .) 

Effect of plan-form shape on the rotational damping-ln-pitch 
derivatives . - A comparison is made in figure 14 between an aspect-ratio- 3 , 
tailless delta wing in combination with a body and an aspect-ratio- 3 , 

52.5° swept highly tapered wing in combination with a body (model 3 ) to 
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determine the effect of plan-form shape on the rotational damping- in -pitch 
derivatives. The data for the delta -wing configuration were obtained from 
reference 5 - 

From figure 14, it appears that there is little effect of plan-form 
shape on the rotational damping-in -pitch derivatives at M = O .915 and 
M = 1.05; however, at M = I.I 6 the delta-wing configuration has better 
damping characteristics than the swept -wing configuration. Reference 1 
also shows that the damping in pitch of the delta-wing plan form is 
superior to the swept-wing plan form at supersonic speeds. 


Drag 

The variation of the drag coefficient with Mach number for each of 
the four models is presented in figure I 5 . The effect of wing sweep may 
be seen by comparing models 3 and 4. Model 4 has lover drag throughout 
the Mach number range than model 3 as would be expected since model 4 has 
a greater wing sweep than model 5. An indication of the effect of plan 
form may be made by comparing models 1 and 2 where it is seen that the 
composite plan form has considerably better drag characteristics than the 
conventional swept plan form even though the wing of the composite plan 
form was somewhat thicker than that of the conventional swept plan form. 


CONCLUSIONS 


From the results of the free-fllght tests of four wing-body combina- 
tions incorporating identical bodies and tested at the same center -of - 
gravity location of I 7 percent mean aerodynamic chord auid consisting of 
an aspect -ratio -4, 45° swept composite wing, an aspect-ratio-4, 45° swept 
highly tapered wing, an aspect-ratio-5, 52 - 5 ° swept highly tapered wing, 
and an aspect -ratio -5, 60 ° swept highly tapered wing, the following 
conclusions may be stated: 

1. All models were statically stable throughout the Mach number range 
investigated (M = O .85 to 1.30). 

2. For the low amplitudes reported in this paper, the total damping 
factor was affected by plan-form shape, whereas the rotational damping- 
in-pltch derivatives of the models were low throughout the Mach number 
range Investigated and differed appreciably between models. Inasmuch as 
the damping-in-pitch derivatives were low, the total damping factor con- 
sisted mainly of the contribution of the slope of the lift curve. 
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3. At M = l.l 6 and for the vrings with aspect ratio 3, the delta- 
wing configuration had better rotational damping-in-pitch characteristics 
than the swept-wing configuration. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Leingley Field, Va., Aiigust 27, 195^- 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF THE MODELS 


Wing: 

Total area, sq ft 

Span, ft 

Aspect ratio 

Mean aerodynamic chord, ft . . 
Sweepback of quarter -chord, deg 

Dihedral, deg 

Taper ratio 

NACA airfoil sections parallel 
to free stream 

Fuselage : 

Length, in 

Fineness ratio 

Miscellaneous : 

Model weight, lb 

Moment of inertia in pitch, 

ly, slug-ft2 

Center -of -gravity position, 

percent M.A.C 

Wing loading, Ib/sq ft .... 
Relative density factor, ^ - 

At M = 0.85 

At M = 1.25 

^•Outboard . 

^Inboard. 


Model 1 

Model 2 

Model 3 

Model 4 

5.76 

4.40 

4.40 

4.40 

4.82 

4.20 

3.64 

3.64 

4 

4 

3 

3 

1.59 

1.21 

1.39 

1.39 

45 

45 

52.5 

60 

0 

0 

0 

0 

ao.715, 'bo. 288 

0.20 

0.20 

0.20 

65A006 

65AOO4 

65AOO4 

65AOO4 

75.00 

75.00 

75.00 

75.00 

10.70 

10.70 

10.70 

10.70 

128.25 

121.00 

125.00 

123.50 

10.96 

10.30 

10.10 

10.50 

17 

17 

17 

17 

22.2 

27.5 

28.4 

28.0 

222 

356 

346 

342 

199 

320 

298 

295 
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(a) Model 1. A = 4; airfoil sections, NACA 65 AOO 6 . 

Figure 1.- General arrangement of models. All dimensions in inches mless 

o*therwise noted. 
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(b) Model 2. A = 4; airfoil sections, NACA 65A004. 
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Figure 1.- Continued. 
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(c) Model 5- A = 5; airfoil sections, NACA G'^kOOk. 
Figure 1.- Continued. 
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(d) Model 4. A = 3; Airfoil sections, NACA 6'^AOOk. 


Figure 1.- Concluded. 
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Figure 2.- Photograph of one of the models. L- 818 I 3 
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Figure 5*- Photograph of one of the models on the launcher. 


L-81915. 
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Figure U.- Longitudinal distribution of cross-sectional area. 
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Time, sec 


(a) Model 1. 

Figure 5*~ Time history of flights showing envelopes of the oscillations 
trim values^ and the times where the pulse rockets were fired. 
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(b) Model 2. 
Figure Continued. 
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(c) Model 5 . 
Figure Continued. 
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1.4 

1.2 

1.0 

.8 



(d) Model 4. 
Figure 5*- Concluded. 
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Figure 6.- Variation of air density with Mach number. 
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Figure 7 .- Variation of velocity with Mach number. 
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Figure 8.- Variation of dynamic pressure with Mach number. 
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Figure 9*- Variation of Reynolds number, based on respective mean aero- 
dynamic chord, with Mach nimnber. 
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(a) Models 1 and 2. 



(b) Models 5 and i+. 

Figure 10.- Variation of the slope of the lift curve with Mach number. 
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(a) Models 1 and 2. 
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(b) Models 5 and 4. 


Figure 11.- Variation of the slope of the pitching -moment curve with Mach 

number. 
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Ca) Models 1 and 2 . 



(b) Models 3 and 4. 

Figure 12 .- Variation of the total damping factor with Mach number. 
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(a) Models 1 and 2. 
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(b) Models 3 and 4. 

Figure 13*- Variation of the rotational damping- in-pitch derivatives with 

Mach number. 
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Figure l4. Effect of plan-form shape on the rotational damping- in-pitch 

derivatives . 
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Fig\jre I 5 .- Variation of the drag coefficient with Mach number. 
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